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Abstract 

We study numerically temperature-shift and field-shift aging protocols on the 3-dimensional (3D) Ising Edwards- 
Anderson (EA) spin-glass (SG) model focusing on respectively the temperature-chaos nature and the stability 
under a static field of the SG phase. The results of the latter strongly support the droplet theory which predicts 
the instability of the SG phase under the field. They are also discussed in relation with the experimental studies. 
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1. Introduction 

More than a decade the dispute on the nature 
of ordering in spin glasses between the mean-field 
theory based on the replica symmetry breaking in 
Parisi's solution of the infinite-ranged model and 
the droplet theory based on a phenomcnological 
scaling argument on the short-ranged model has re- 
mained unsettled [1] . Currently, the following two 
explicit predictions by the droplet theory [2,3] have 
been extensively studied. 1) The equilibrium SG 
states at two temperatures with difference AT are 
uncorrclatcd in length scales larger than the over- 
lap length Lat, i-e., the so-called temperature(T)- 
chaos nature of the SG phase, and 2) in the equi- 
librium and thermodynamic limits the SG phase is 
broken by a static field h of infinitesimal strength, 
thereby introduced is the crossover length Lh- The 
latter separates characteristics of droplet excita- 
tions by their size L, such that they are dominated 



by the Zeeman energy for L > Lh and by the SG 
free energy gap for L < Lh- 

One may imagine that the existence of such 
a characteristic length Lat (Lh) could be easily 
checked by the temperature(T)-shift (fields- 
shift) aging protocol, where the temperature 
(field) is changed once on a way to equilibrium. 
However, this is rather hard because of the ex- 
tremely slow spin dynamics in the SG state: within 
a time- window of relative magnitude of 10 4 ~ 5 , 
which is the case not only in simulations but 
also in ordinary experiments, the SG correlation 
length, or mean domain size, R(t), is considered 
to grow only by a factor two or even less as will 
be discussed below. Therefore, no simulation nor 
experiment on such an aging protocol alone has 
not yet succeeded to settle the problem mentioned 
above. A possible strategy to overcome this diffi- 
culty is to examine a set of T- (h-) shift protocols 
with Lat (Lh) systematically changed relatively 
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to R(t) within the available time-window, to look 
for a possible scaling behavior in the length scales 
involved, and to check whether the behavior is 
compatible with the droplet theory or not [4] . 

In the present work we numerically study ag- 
ing dynamics in the 3D Ising EA SG model with 
Gaussian nearest-neighbor interactions with zero 
mean and variance J. The latter is used as the 
unit of energy and temperature. In this unit the 
SG transition temperature is estimated as T c ~ 
0.95 [5,6]. Along the above-mentioned strategy we 
have examined a number of T- and /i-shift pro- 
cesses with various AT and h as well as the (wait- 
ing) time, i w , of the shift. For further details of the 
model and the Monte Carlo (MC) method of sim- 
ulation the reader may refer to our previous works 
[7,8,9,10,11]. 

In the next section we investigate the T-chaos 
nature through T-shift processes, and in §3 the sta- 
bility of the SG phase under a field through /i-shift 
processes. The last section is devoted to discus- 
sions. 



2. T-shift aging processes 

Following the idea of 'twin-experiments' in [4], 
we examine here negative (positive) T-shift aging 
processes. A system is instantaneously quenched 
at time t = from T = oo to T (T 2 ) below T c , 
the temperature is decreased (increased) to T 2 (T) 
at t = (£P), and then let the system to age at 
T 2 (T). Throughout the process we measure the 
imaginary part of the ac susceptibility of frequency 
u), x"{u;t) [8]. Its typical behavior in a negative 
T-shift process is shown in Fig. 1. By the T-shift 
at t = t™, x"( w ; t) rapidly decreases, and it merges 
to the T 2 -isothermal curve from below. If, how- 
ever, the branch of at t > is shifted 
to the right by a proper amount, t s h, as indicated 
by the open circles in the figure, it merges to the 
T2-isothermal curve from above at t = t mr . The 
inset of the figure explains how to determine t s h- 
when an amount of the shift is smaller than t s h, 
the shifted branch crosses with the T^-isothermal 
curve (triangles in the inset), while a larger shift 
delays the merging significantly (squares). Here we 
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Fig. 1. x" (lo\ t) in a negative T-shift process with = 512, 
Ti = 0.7, T 2 = 0.4, and t w = 2tt/uj = 64 (solid circles). The 
two continuous curves represent x"(u>;i) in the isothermal 
aging at T = Ti and T2. The open circles are the properly 
shifted branch of x"(u;;t) at t > t° . Shown in the inset 
are how the shifted branches merge to the T2-isothermal 
curve (see the text). 
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Fig. 2. (t^,) vs. tJJ, (tcr) for the negative (positive) T-shift 
processes between three pairs of temperatures. The lines 
represent the expected behavior from Eq.(l) combined with 
Eq.(3) as explained in the text. 

call the time required for x"( w ; t) to merge to the 
isothermal curve in this sense the crossover time 
and denote it as i" r , i.e., i" r = £ mr - (t™ + t sri ). 
[This characteristic time is called the effective age 
in [4] and the effective waiting time in [11].] For a 
positive T-shift process we can similarly define the 
crossover time i[? r [11]. A set of results for t w and 
ten thus obtained are shown in Fig. 2. 

After the first quench of a negative T-shift pro- 
cess, the SG ordering, or the mean domain size 
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R(t), grows up to R(t w ) at t — t w . By the T-shift 
R(t) does not exhibit a discontinuous change, but 
continues to grow now with the growth rate at T2 
[9]. However, it takes time t CI for x"(o>; t) to merge 
to its T2-isothermal curve. We suppose that in this 
time range, which we call the transient regime, sub- 
domains of a mean size L(t — t w ) in local equilib- 
rium of T2 grows within each domains of a mean 
size R(t w ) in local equilibrium of T\. If R(t w ) is 
enough smaller than the overlap length Lat in- 
troduced in §1, we naively expect that L(t — t w ) 
catches up R(t w ) at t — t w ~ i cr , or, 

L(t CI ) = R T2 {t CI ) = R Tl {t w ), (1) 

to hold, where Rr{t) is the growth law of SG do- 
mains in the isothermal aging at T. By the first 
equality above we imply that the growth mecha- 
nisms of L(t~t w ) at T and Rr(t) are common, i.e., 
they are due to thermal activation processes. The 
same equation as above but with T\ and T2 inter- 
changed is expected to hold for the corresponding 
positive T-shift process. This expectation is in fact 
confirmed by the result in Fig. 2 that the data t" r 
vs. t^ and t^ vs. t£ r lie on a common curve for a 
small AT = Ti - T 2 (= 0.1). We call the inter- 
pretation so far described the cumulative memory 
scenario for T-shift processes with small |AT|. 

For later discussions let us here mention explicit 
expressions for the growth law of Rx(t). Within 
the accuracy of simulations, the numerical results 
of R T (t) [7,12,13] are well fitted cither to the log- 
arithmic law proposed in the droplet theory [14] 

Rr(t)/lo ~ l^Ht/to)]^, (2) 

or to a power law 

i? T (i)A>^Ao) 1/z(T) . (3) 

Above, Iq and to are the characteristic length and 
time scales, and A g and ip the characteristic scale 
of energy barrier of droplet excitations and the as- 
sociated exponent, respectively. In Eq.(3) we set 
?o = 1 lattice distance and t = 1 MC step per 
spin, and the exponent l/z(T) is given by [7] 

l/z(T) ~ 6T/T C , (4) 

with b ~ 0.16 in the range 0.7 > T/T c > 0.4. 



For a large AT (= 0.3), the sets of data of neg- 
ative and positive T-shift processes in Fig. 2 dif- 
fer significantly from each other. Also the second 
equality of Eq.(l) is violated for the negative T- 
shift process shown in Fig. 1. We tentatively at- 
tribute this deviation from the cumulative mem- 
ory scenario to the temperature-chaos nature of 
the equilibrium SG phase described in §1. If Lat is 
smaller than Rx t (£ w ), then in the time range after 
the T-shift in which L(t — t w ) > Lat is satisfied, 
aging dynamics looks as if the system is already 
in the T2-isothermal aging state. This is because, 
by definition of Tat, the longer range order than 
Tat developed at T = T\ is irrelevant to the SG 
order at T = T2. Therefore t cr is given by the con- 
dition L(t cr ) ~ Tat, and is independent of t w in 
this case. 

However, t w in the present simulation is not large 
enough to reproduce such a drastic phenomenon. 
We consider that the negative T-shift process ob- 
served here is in the weak chaos regime proposed 
by Jonsson et al [15]. We also note that the second 
equality of Eq.(l) is little violated in the positive 
T-shift with AT — 0.3. This asymmetric appear- 
ance of the chaos effect contradicts the prediction 
of the droplet theory [2] as well as the recent ex- 
periment on the Heisenberg spin glass AgMn [4] . 



3. h-shift aging processes 

The /i-shift processes we mainly discuss here are 
as follows. After instantaneous quench to T, the 
system is aged under a zero field until t = t w , 
when a field h is switched on. We then observe the 
induced (zero-field-cooled) magnetization M (= 
Mzfc) an d x"( w ; t) in the whole process. We study 
them exclusively at T = 0.6 where the aging dy- 
namics is considered not affected by the critical 
dynamics close to T c . 

It is well established that, at least for a small 
h, the logarithmic growth rate of M defined by 
S(t') = d{M/h)/d\nt' with t' = t - t w exhibits a 
peak at t' ~ t w . More generally, the peak position 
of S(t') plays the same role as t cr discussed in the 
previous section. It gives us a characteristic time 
scale of crossover from the isothermal aging state 
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Fig. 3. S(f) with t w = 4096 plotted vs Inf. The data arc 
for h =0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.75, 1.0 and 2.0 from 
right to left. 
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Fig. 4. Plot t cr vs. t w . The data points are for h =0.1, 0.2, 
0.3, 0.4, 0.5, 0.6, 0.75 and 1.0 from top to bottom, and the 
line represents i cr = t w . 

under h = to that under a finite h, and so is 
denoted also by t CT . In Fig. 3 we show S(t') with 
t w — 4096 for various values of h. For a small h (= 

0. 1), as mentioned just above, i cr ~ t w is ascer- 
tained. When h increases, however, t cr significantly 
decreases. In Fig. 4 the resultant t cr for various h 
and t w are presented. 

Now we convert the relations between t cr and 
t w into those between the characteristic length 
scales R^f — Rr(t cr , h) and R T (t w ). The effective 
crossover length R^f is the mean size of subdo- 
mains grown under h in the period of i cr and has 
to be estimated by the growth law under h. Since, 
however, it is quite time-consuming to simulate 
the latter due to the presence of the paramagnetic 
component, we approximate it by RT(t C r,h = 0), 

1. e., by Eq.(3) in the present work. 
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Fig. 5. R^/L h vs R T (t w )/L h . The data at small x are 
fitted to y = x — ca; 1 + 1 /' 5 with c = 0.15, which is the 
formula proposed in the weak chaos scenario in [15]. 

The obtained sets of R^ and i?r(iw) for differ- 
ent h are normalized by the crossover length Lh 
introduced in §1. It is given by [2] 

-s 



L h /i ~ (h/ry 



(5) 



with S = (| — O)^ 1 , where T is the unit of free- 
energy gap. In evaluating explicitly Lh by Eq.(5) 
we simply put T unity and 5 = 0.77 with d = 3 and 
9 = 0.2 [7,16]. The resultant scaling plot of R c c ?/L h 
vs. RT(t w )/Lh is presented in Fig. 5. One sees in 
the figure that most of the data tend to lie on a 
single curve. If we take into account the fact that we 
have not adjusted the parameters involved at all, 
the scaling behavior obtained is quite satisfactory. 

We interpret the scaling behavior shown in Fig. 5 
as dynamical crossover from the SG aging state 
to the paramagnetic state. Actually, in the range 
i«l, where Rr(t) measured within the time win- 
dow of our simulation is much smaller than Lh, 
various properties which are usually regarded as 
those of the SG state are observed. Among them 
are the so-called wTj-scaling of the ac susceptibil- 
ity in the isothermal aging [1] and the sum rule 
M ZFC (t) + Mtr(*) = Mfc(*) [17], where M FC is 
the induced magnetization observed when h is ap- 
plied just after the quench to T and M TR is the 
thermoremanent magnetization observed after h in 
the field-cooled process is switched off at t = t w . 
In the range x ^ 1 corresponding to a large field of 
h ^ 1.0, on the contrary, the complete saturation 
of Mzfc and Mfc to a common value is observed 
within the time window of our simulation. This is 
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the expected property of the paramagnetic state. 
Combined with these observations which will be 
discussed elsewhere, we conclude that the result of 
Fig. 5 is a strong support of the droplet picture for 
the /i-shift aging process and so the instability of 
SG state in the equilibrium and thermodynamic 
limits. 

Lastly let us try to extend the scaling behavior 
obtained above by our simulation on the Ising EA 
model to the one in real Ising spin glasses such as 
Feo.sMno.sTiOs, focusing on the crossover point 
defined by the condition x = i?r(iw)/ Lh = 1- Here 
we use the power growth law of Eq. (3) with to = 
10~ 12 s. Then, with the same approximations, such 
as Rr(t, h) ~ Rt{T, h = 0), as those already in- 
troduced, the condition is reduced to 

h~ s c (t/t )- bT ^. (6) 

Also using b — 0.16 and S = 0.77 as in the simula- 
tion, we obtain, for example, h ~ 0.6T for T/T c = 
0.6 and t = 10 2 s. A set of these figures is com- 
pared with h ~ 1.3 T, Tat/T c ~ 0.6 and t mc ~ 
10 2 s, which are picked up from the paper by Katori 
and Ito [18], who determined the critical temper- 
ature Tat under h of Feo.sMno.sTiOs as the one 
where the Mzfc starts to deviate from Mpo In a 
certain range of Tat, they found the relation h oc 
(1 — Tat/Tc)" with a ~ 1.5 which is the expected 
behavior of the de Almaid-Thouless phase tran- 
sition [20] in the infinite-ranged SG model. How- 
ever, the rather good coincidence of the above two 
sets of vaules, combined with the dependence of 
h on T from Eq.(6) with a fixed t equal to the 
measuring time t mc in their temperature sweeping 
experiment (not shown), suggests that Tat they 
measured could be interpreted as the dynamical 
crossover point with t mc ~ t cr . 

Another experiment on Feo.sMno.sTiOs under 
large static field h was done by Mattsson et al. 
[19] , who measured the ac suceptibility. They found 
two regimes in its behavior which are separated 
by the the freezing temperature Tf (w, h). Although 
what they measued is droplet excitations almost 
in equilibrium but not domains in aging, we may 
use Eq.(6) with t = t w = 2tt/lj, i.e., Tf is the tem- 
prature at which the mean size of droplets which 
can respond to the ac field of period t u is equal 
to L h . Then Eq.(6) yields, for example, h ~ 3.0 T 



with t u = 10~ 2 s and T f /T c = 0.5, while the cor- 
responding experimental value is h — 2.2 T. 

Although we do not want to insist strongly the 
semi-quantitative agreement of the experimental 
results with a simple extension of the simulational 
results, we may claim from the above argument 
that t CT for field-shift processes under ft of a few 
Oe is beyond the astronomical time. 



4. Discussions 

On the interpretation of aging phenomena in 
spin glasses by means of the droplet picture, the 
length scales of certain domains and droplets play 
a central role. In the above arguments, they are 
Tat, Lh, Rr(t w ), and L(t — t w ). In the numerical 
study so far carried out [7,12,13], Rr(t) is replaced 
by the coherence length of the real replica overlap 
function which has been explicitly evaluated. Since 
the overlap function is obtained by averaging over 
sites in each sample and then over samples, de- 
tailed properties of each domains, such as whether 
they are fractal or compact, are averaged out, and 
so they are beyond the scope of the present simula- 
tion. Neither the nature of subdomains has been di- 
rectly pursued by simulation yet. But its introduc- 
tion with the subdomains-within-domain scenario 
[9,11] helps us very much to interpret the observed 
phenomena in the transient regime after the T- or 
/i-shift as well as the memory effect observed in the 
T- or ft-cycling processes [1]. For example, even 
in the cummulative regime of Fig. 5 with R c c f ~ 
-Rt(^w), rejuvenation-like behavior (or a jump-up) 
in x"(<^;£) has been observed both in experiment 
[21] and the present simulation (not shown). 

Another fundamental ingredient of the analysis 
based on the droplet picture is the growth law of 
Rx(t), which is frequently used as the conversion 
formula from the time scale to the length scale. 
However it is hard to accurately determine it due 
to the following circumstances. Let us consider the 
ratio r = i?T (imax)/.RT(£min), where i max (t min ) 
is the maximum (minimum) time scale of observa- 
tion. The results of our simulation with i m i n ~ 10 
and i max ~ 10 5 yield, say at T/T c = 0.6, r ~ 2.3. If 
Eqs.(3) and (4) with t a = 10~ 12 s and b = 0.16 are 
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extended to the time scale of ordinary experiments, 
i.e., i m in ~ 10 s and i max ~ 10 5 s, the identical r 
is obtained. When the logarithmic law of Eq.(2) is 
applied to the experimental data, r becomes even 
smaller [4]. We may say that, both in experiments 
and simulations, we measure a very limited por- 
tion of an equilibration process, at least the length 
scale of SG ordering is concerned. This is the rea- 
son of difficulty in determining the growth law of 
Rr{t) accurately. 

In order to investigate the stability of the equi- 
librium SG state against perturbations such as 
the temperature change AT and a static field h, 
thereby overcoming the above-mentioned diffi- 
culty, we have carried out a set of shift- aging pro- 
tocols where the perturbation of various strengths, 
P (AT or h), is applied after the waiting time t w 
which are also systematically changed. For each 
process specified by (P, t w ), we have extracted the 
crossover time t cr between the state aging to an 
unperturved SG state and the one to the perturbed 
SG, or non- SG state. The length scales R^(P; t„) 
associated with t CT are then found to be scaled 
nicely by the characteristic length scale associated 
with P. This strategy, first carried out in [4], is 
expected to be most efficient to judge the stability 
of the equilibrium SG state by actual experiments. 

Irrespectively of the strategy mentioned above, 
experiments to extract Rr{t) with referring to 
another length scale independently of the growth 
law itself are certainly very much required. To our 
knowledge, however, there has been only two such 
experiments which use a length scale associated 
with the static field effect. One is due to Mattsson 
et al [19] already mentioned above and the other 
due to Joh et al [22]. The latter result prefers the 
power law of Eq.(3) with the value of b in l/z(T) 
nearly equal to our numerical results. From this 
point of view, an experiment on SG fine particle 
systems may be of interest. 

To conclude we have presented and discussed the 
results of our simulation on the T- and /i-shift aging 
processes in the 3D Ising EA SG model. In the T- 
shift process only a precursor of the temperature- 
chaos effect has been observed, while the results 
on the /i-shift process strongly support the droplet 
picture that the SG state under a finite field is 
unstable in the equilibrium and thermodynamic 



limits, though it takes astronomic time for the SG 
state to be broken if h is small. 

I would like to sincerely thank K. Hukushima, 
P. Jonsson, H. Yoshino, and E. Vincent for their 
fruitful discussions. The present simulation has 
been performed by using the computor facilities 
at Supercomputer Center, Institute of Solid State 
Physics, the University of Tokyo. 
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